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PREFACE 


This  RpOft  IS  a  summary  of  the  fmal  portion  of  the  woik  performed  on  this  contract  which  deals  with 
modeling  die  BSD-caused  igubon  of  solid  propellants.  A  summary  of  the  pnor  woric  was  already  submitted  in  an 
interim  report  that  preceded  an  ARO-sponsoted  woricsbop,  held  in  NashviUe»  Tennessee  on  Apnl  17  and  18. 1989, 
devoted  primarily  to  the  question  of  igoitioa  following  an  BSD  breakdown  event  That  mtenm  report  dealt  only 
poqifaetiliy  with  ignitioa 

However,  there  was  always  an  intent  to  consider  ignmon  in  the  final  year  of  the  contract,  so  this 
docuraeot  adds  materiftl  00  that  subject  as  well  as  contaimog  several  new  items  of  research  not  reported  m  a  thesis, 
that  is  the  other  mam  ouqmt  of  this  contract 

The  authors  wi^  to  tfmnk  ARO's  Projea  momtor,  Dr  David  Mann,  for  taking  a  strong  mterest  m  the 
subject  matter  of  this  project 
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CHAPTER  1 

EXECUTIVE  SUMMARY 


This  docomem  is  the  ftoal  rqxnt  for  t  three  year  study  performed  at  Electromagnetic  Applications,  Inc 
(EMA)  for  the  U^.  Anny  Research  Office  (ARO).  The  subject  matter  is  an  analytical  treatment  of  microstructural 
tipects  of  Electroshuic  Discharge  (ESD)  ignition  of  solid  rocket  motors.  This  report,  and  a  companion  document 
(Dr.  Marie  Cyore's  thesis)  [1],  are  an  expansion  of  an  intenm  report  covering  similar  matenal,  released  in  mid  1989 
[2]. 


ESD  hazards  of  the  type  described  herein  are  now  well  known  to  the  solid  propellant  industry  [2],  although 
this  was  less  so  the  tune  of  our  tuiginal  juoposai  to  ARO  m  1986.  Since  1985,  two  accidents  involving  solid 
rocket  motors  have  die  loss  of  eight  bves.  These  are  now  accepted  to  have  been  imtiated  by  ESD,  causing 
ESD  to  have  become  the  leadmg  safety  hazard  wuhin  the  sohd  propellant  mdusoy; 

This  report  describes  features  of  solid  propellants  that  make  them  particularly  suscepu1)le  to  ESD  Primary 
emphasis  has  been  placed  on  nucrostnictural  modeling  to  better  understand  how  to  mitigate  these  hazards  by 
undosUDidiiig  die  details  of  die  hazard. 

Attempts  have  been  made  m  this  and  the  companion  doctoral  thesis  [1]  to  explain  (among  other  things) 
effects  of  tune  (or  frequency),  propellant  formulabon,  temperature,  pressure,  and  relauve  humidity  This  report 
mqilaiiis,  but  is  not  intended  as  a  final  definitive  explanation  of,  mudi  of  the  experimental  phenomena  that  is 
sommanzed  herein,  in  summary,  we  can  make  the  following  general  statements  that  are  supported  by  the  (primarily) 
micioscc^ic  analyses  within  this  ARO  contract 

1.  As  wmk  has  been  performed  on  the  electrostatic  breal^own  and  ignition  of  sohd  propellants  m 
response  to  the  two  accidents,  it  has  become  clear  that  there  ate  good  macroscopic  and  microscopic 
reasons  for  ESD  to  be  a  senous  hazard  for  sohd  propellants.  Based  on  both  static  and  transient 
macroscopic  sunulations  by  EMA  (not  the  subject  of  this  report),  the  propellant  fields  from  realisuc 
charge  densioes  can  sometimes  exceed  those  known  to  cause  breakdown  in  small  expenmental 
sanqiles.  Fnxn  a  modeling  standpomt,  the  analysis  of  the  fi^  distnbution  is  not  simple,  but  good 
accuracy  is  possible,  and  there  are  few  doubts  that  hazardous  situations  can  occur 

2.  From  the  imeroscopic  viewpouit,  the  general  features  are  also  understood  Microscopic  modeling  is 
able  to  explam,  prunarily  through  the  stausucal  details  of  the  proximity  of  the  alummum  parucles, 
many  of  the  expenmental  breakdown  variations.  However,  the  details  of  the  behavior  of  propellant 
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bmdcdown  fields,  condocuvity  and  pemituvny  as  a  funcnon  of  tune  (or  frequency)  and  temperature 
are  only  incompletely  undctstood. 

3.  Among  the  vanous  mflcrnmodeling  techniques,  the  fuute  element  technique  has  shown  the  most 
stgnificant  advantages  in  terms  of  accuracy,  case  of  mteracuon  dunng  problem  set-up,  case  of 
obtaining  solution,  quality  of  display  and  mteipretanon  of  results,  short  computing  time,  and  total 
solncioo  time.  However,  the  fmite  element  method  has  not  been  of  significant  value  m  studying 
random  medu  (exc^t  as  u  occasionally  is  used  as  in  the  above  ways),  since  it  is  not  usually  used  for 
multiple  runs  with  vasyrng  parametersfsnch  as  m  sequential  brealcdown  discussed  in  [Ref  1]  For 
this  study  of  random  media,  resort  is  largely  made  to  boundary  element  methods  [1], 

Among  the  other  oooclusnns,  mostly  discussed  m  detail  m  [Ref.  1]: 

Bredaloini  fields  (at  least  in  the  cases  studied  to  date)  are  Imearly  related  to  average  mterpamcle  spacmgs  m 

both  highly  Structured  and  random  media.  The  relauonship  is  simpler  than  ongmally  expected. 

Formulas  have  been  obtained  for  interparticle  fields  between  spheres  with  accqnable  accuracy  and  speed  of 

sdatioa. 

New  insights  have  been  gamed  on  the  effea  of  differential  conducuvity  in  explaimng  both  field 

enhancement  and  resistivity  vs.  tune  experiments. 

Much  of  the  materiai  that  would  normally  appeared  m  this  final  report  have  already  appeared  m  Reference 
[2],  and  vutnally  none  will  therefore  be  rqieated  here.  A  brief  overview  was  given  m  Secuon  2  of  [2]  of  the  four 
major  parts  of  the  breakdown  phenomenon.  This  was  followed  m  us  Section  3  by  a  descnpuon  of  a  number  of 
modeling  Kvhnimiw!  and  results  that  are  prunarily  useful  m  macroscopic  modeluig,  but  with  emphasis  on  then-  use 
m  microscopic  mmteling  Sections  4  and  5  dealt  with  major  poruons  of  Dr  Cyure's  thesis,  which  now  appear  as  a 
separate  componeut  of  this  final  report  [1]  Brief  comments  were  made  m  Section  6  of  a  numba  of  topics  of 
imptutance  to  undeislandmg  the  physical  aqiecis  of  breakdown  not  covered  m  the  thesis;  this  woik  was  performed  by 
ofiier  mvestigatais  rvmkmg  on  the  projecL  The  final  secuon  of  the  mierim  report  [2]  gave  some  of  the  team's  views 
on  items  diat  were  covered  m  the  ARO-sponsored  workslx^  [3]  for  which  the  mterun  report  was  prepared  That 
secfioo  has  been  reiamed  and  expanded  as  pan  of  this  document 

A  major  poruon  of  the  team  effort  was  the  support  of  Mr  (now  Dr.)  Mafic  Gyure  m  obtaming  h,s  doctoral 
degree  (m  Physics)  at  the  Umversity  of  Colorado  That  work  was  under  the  direcuon  of  Dr.  Paul  Beale,  who  has 
also  authored  some  sections  m  the  present  volume.  Mark  Gyure's  doctoral  work,  in  its  entirety,  is  contained  m  the 
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Kpnte  lepott  [1].  He  has  also  contributed  a  short  section  on  energy  considerations  related  to  his  thesis  work  in 
this  final  report 
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CHAPTER  2 

PARTICLE  SEPARATIONS  AND  EXPERIMENTAL  VERIFICATION 
(Prepared  by  Dr.  Paul  Beale) 


One  of  the  piunacy  lesults^  of  the  elecoostatic  bicakdown  modeling  is  the  simple  connection  between  the 
didectnc  breakdown  streugth  of  a  metal-loaded  dielectnc  and  the  average  nearest-neighbor  surface-to-surface  distance 
between  die  metallic  particles  m  the  composite. 


Eb  =  Eb'>^;i 


(2.1) 


The  breakdown  field.  Et>.  is  directly  proporuonal  to  the  bre^own  field  of  the  binder.  The  proportionality  constant 
dqiends  on,  d,  the  average  nearest  neighbor  distance  between  metal  particles  and.  D,  the  diameter  of  the  metal 
particles.  Unfortunately,  die  distance  between  particles  m  a  random  composite  is  difficult  to  measure.  Instead  we 
d^ennined  this  value  by  ifipeaiing  to  known  mathematical  pn^iemes  of  the  correlauons  between  hard  sphencal 

ia  die  mechanics  of  single  fluids. 

Let  p  ss  N/V  be  die  number  densuy  of  pansies  in  a  volume  V.  The  pair  ctMTclation  funcuon  g(r)  is  defined 
by  4  It  p  g(r)  dr;  tt  is  the  probability  that  the  center  of  a  pamcle  is  located  at  a  distance  between  r  and  r-Klr  from 
the  center  of  a  chosen  particle.  The  average  nearest  neighbor  distance  can  be  determined  fiom  the  dismbuuon  funcuon 
for  particle  sqiaratioQS.  Let  p(r)  dr  be  the  probability  that  the  closest  neighbor  of  a  given  pamcle  is  at  a  distance 
between  r  and  r^dr  fiocn  a  given  particle.  It  is  not  too  difiicult  to  show  that 


p(r)  =  4itpr^g(r) 


(2.2) 


From  this  distriboaou  function,  the  detennmation  of  d  is  easy. 


d=Jrp(r)dr-D  (2.3) 

0 

Now.  the  goal  of  the  nucrostructuial  theory  is  to  detemune  the  pair  correlauon  funcuon  g{r)  The  added 
difficjlly  is  the  situation  we  have  with  propellants  is  that  the  composue  is  composed  of  different  types  and  sizes  of 
pamcles.  Let  Ni  be  the  number  of  particles  of  type  i  with  diameter  D,.  The  number  densiues  of  the  different 
components  are  given  by  pi=Ni/V.  There  wdl  be  several  different  pair  cotrelabon  funcuons  correspondmg  to  the 
coirelanons  between  paindes  of  type  i  and  pamcles  of  type  j.  Let  p,  pj  gij(ri-r2)  d^i  d^r2  be  the  probability  that  a 
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pvude  if  type  i  is  located  inside  volume  element  d^ri  located  at  posiuon  ri  and  a  parucie  of  type  j  is  located  inside 
volume  element  d^  located  at  position  t2.  The  nearest  neighbor  probability  distnbuuon  funcuons  are  pjjO)  dr  is  the 
probability  diat  t  panicle  ^pe  i  has  its  nearest  neighbcM*  of  type  j  at  a  distance  between  r  and  r-i-dr  This  is  gotten 
foun  the  pair  COTcela&oo  functions  gij(r)  by  the  relauon 


pij(r)a4npjgij(r) 


(2-4) 


The  goal  of  oar  miciostnicntral  theory  is  to  dcien>..i'^  ‘  ~ir  fane  and  tro.n  those  to  deiermine  the 

nearest  neighbor  ^pacing  between  metal  paraclc:»  The  n'tOde:  wdi  be  rsstmed  is  t^lat  the  pa.'tK'ies  arc  hard 

^pboes  of  diameter  Di  widi  number  denaues 

The  calailatmn  of  the  pair  conelatton  funcuons  proceed:,  the  t  oltzmarn  assunipur  tiiat  any 
configuration  of  particles  vduch  is  possible  (no  hard  spheres  o*  erlapj  has  eq  ial  a  pru^r  p.*obabtiii\  T'nce  th>s 
also  the  basis  of  all  of  statistical  mechanics,  we  wul  ma^  use  of  the  stausucal  mechai..ca  of  a.nri'  uquids  to 
deteniuDe  die  pair  correlation  funcuons.^  The  particular  approach  used  is  based  on  uc  properues  of  the  direct 
correlation  fimccioas  qj(r)  which  are  defined  by  the  Omstem'Zernike  equation.^ 

g,j(r)  - 1  =  c,j(r)  +  Pit  /cik(r') !  gkj(r-r’)-l  ]  dr  (25) 

1; 


This  10  Itself  contains  no  new  infonnauon.  The  calculauon  of  the  pair  correlation  funcuons  depends  on  being  able  to 
close  diese  equatioos  with  a  means  of  specifying  another  constraint  on  the  direct  correlauon  funcuons.  The  constraint 
we  use  is  the  Perkus-Yevick  equauon*^'^  ^  which  for  hard  spheres  states  that  the  direct  correlauon  funcuon  Cyfr)  is 
zsm  fbrallr>(D|  +  Ilj)/2.  This  constramt  allows  the  exact  calculauon  of  all  of  the  direct  ccorelauon  funcuons  and 
pair  corroiaDon  fimetioas  for  up  to  two  different  ^here  sizes  This  approxunauon  is  well-tested  m  the  stausucal 
mechames  of  simple  liquids  and  is  known  to  be  quite  accurate  even  at  very  high  panicle  densiues  near  the  phase 
transition  toa  solid  crystalline  phase  ^ 

Lebowitz  gives  an  exact  expression  for  the  direct  correlauon  funcuons  ^  ^  They  are  of  the  form 
r  ai  +  bj  r  +  d  r  <  Di 

-cu(r)=i  (2  6a) 

I  0  r>Di 
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I' 


32  +  b2  r  +  d  r  <  D2 


>e22(r)  = 


-C12{r)«^ 


*1  + 


0  r  >  D2 

31  r<l 

b  +4  1  J  3^  +  d  3^ 


l<r<Di2 
r>Di2 

where  3=r-l  3nd  the  coeffidents  si .  bi,  b  3nd  d  3re  given  in  Lebowitz's  paper  ^  ^ 

The  Fourier  transforms  of  the  dntct  conekition  {imcbons  3ie  defined  by 

«0 

Cij(k) = Y  ^  ^  ^  ^ 

0 

and  ttke  Foaner  transforms  of  the  pair  coodaDon  funcQon  arc  defined  by 


0 


In  Fourier  transform  space  the  Omstem  Zenuke  equation  has  the  form 


(2  6b) 


(2.6c) 


(27) 


hijCk)  =  c,j(k)  + 1  p„  C,a(k)  hnj(k) 


this  equation  can  be  solved  for  the  Founer  transforms  of  the  pair  coirelauon  funcuons  hij(k)  Finally  the 
determination  of  the  pen-  coitehmoa  functions  is  accomplished  by  inverse  Founer  transforming 


(2  8) 


(2.9) 


Eij(f)  =  1  +  k  h,j(Ii)  sin{  kr  )  dk 


(210) 


Note  that  this  step  cannot  be  done  analyncaUy  (at  least  not  easily)  so  this  final  step  must  be  done  numcncally  using 
a  fast  Fourier  transfonn  routme.  Once  this  is  accomplished  equations  (3)  and  (4)  can  be  evaluated  numcncally  to 
determine  the  nearest  neighbor  spacmgs  in  tne  composite 
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Our  Uieoietical  claim  in  equaaon  (1)  that  the  breakdown  field  scales  simpiy  with  the  nearest-neighbor 
placing  can  be  tested  with  expenmental  data.  Table  2.1  shows  the  mixture  proportions  and  breakdown  field  for  a 
senes  of  men  propellants  mixed  by  Gibson  at  Edwards  AFB  and  tested  by  Covino  and  Hudson  at  China  Lake 
NWC.^^  In  the  samples,  the  ammomum  perchlorate  was  replaced  by  NaCl.  From  the  mixture  proportions  the  pair 
sozrelation  funcdtms  and  Dearest-neighbor  ^lacmgs  were  determined  using  the  above  procedum  The  NaCl  particles 
were  assumed  to  be  100  nmeR  larger  tnan  the  A1  particles.  Figures  2  1  -  2.3  shows  the  pai  orrelation  functions  for 
sample  KJ-IS  and  Figure  2.4  shows  the  Al-Al  nearest  neighbor  dtsmbuDon  function  The  calculated  nearest  neighbor 
qiacuig  fOT  this  sample  was  0.04S  Dj .  Figare  2.5  shows  the  experimental  breakdown  field  plotted  vs  d/CD^j  +  d) 
Equatum  (2.1)  predicts  that  the  should  fall  on  a  straight  line  with  slope  intersect  the  ongin  As  is 

evident,  the  data  is  well  qiproximated  by «.  straight  line  (whidi  mteisects  the  ongm  as  predicted)  and  that  the  slope  of 
19MV/m  IS  m  goo&  agreement  with  the  subsequently  determined  expenmental  bmdei  breakdown  field  of  23KfV/m 
Ihis  mtcelieoi  agreement  between  theory  and  &.tpenment  is  vtry  encouragmg  and  constitutes  a  quanutauve  theoretical 
approach  for  A*i<^inrr>g  the  dielecmo  brea’tdown  strength  of  random  metal-baded  dielectric  composites 


Table  2.1 

Mixture  Proportions  and  Breakdown  Fields  for  Inert  Propellants 


Sample 

Volume  % 

3  micron  A1 

Volume  % 
NaCl 

Breakdown 

Field  fMV/niy.... 

d/D 

KI-17 

3.0 

55.7 

4.7 

0.327 

KJ-14 

1Z7 

47.5 

165 

0.087 

KI-15 

22.4 

35.9 

0  75 

0  045 

KJ-16 

334 

25  0 

045 

0  025 

KJ-lg 

54.7 

0.0 

0.25 

0.0126 

I 

The  expenmental  breakdown  field  plotted  versus  the  theoretically  determined  average  nearest  neighbor 
n  dustance.  The  data  is  shown  in  Table  2.1.  The  slope  of  the  line  is  18  MV/m  which  is  the  theoretical  prediction  of  the 

*  breakdown  &Id  die  binder.  The  subsequendy  determined  breakdown  field  for  the  binder  was  23  MV/m 

I 

I 

I 

I 
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Figure  2.2  The  NaCl'NaCI  Pair  Correlation  Function  for  KJ-15 


g12(r) 


«  Breakdown  Field  (MV/m) 


0  0.05  0.1  0  15  0.2  0.25  0.3 

d/(d+D) 


Figure  2.5  Test  of  Theory  on  China  Lake  Data 
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CHAPTER  3 

ENERGY  FOUND  IN  MlCROSCOHC  MODELING 
(Prepared  by  Dr.  Mark  Gyure,  as  an  Extension  of  Work  in  his  Thesis  [1]) 


Hk  wolk  of  Reference  [1]  is  pnmanly  addressing  the  problem  of  breakdown  in  randomly  inhomogeneous 
maerials  and  rot  the  problem  of  igniticm  of  acmai  materials  such  as  solid  rocket  propellants.  The  breakdown 
simnlatioos  described  in  this  woik  do  provide  enough  data,  however,  to  estimate  the  energy  released  in  local 
breakdown  events  and  this  information  should  be  useful,  if  not  essential,  to  anyone  interested  in  modeling  ignition 
processes  m  these  mataials. 

The  ignition  problem  of  interest  to  the  aerospace  community  can  be  stated  in  the  following  way  If  a  solid 
rocket  propeDant  expen«ices  dielectric  breakdown,  what  is  the  probability  that  the  propellant  will  also  ignite  and 
bant?  Or  put  another  way.  is  electrical  breakdown  a  suHicient  condition  for  igmuon  of  a  given  propellant  or  can 
bceakdown  occur  widiout  ignition?  The  problem  of  tgniQon  of  soha  rocket  propellants  due  to  dielectric  breakdown  is 
a  vtty  difiicaU  oat  foi'  several  reasons  that  will  not  be  discussed  in  detail  here.  Bnefly.  the  problems  encountered  are 
the  inhoaiogeaeity  of  the  maienal  and  the  difficult  of  modeling  a  system  of  sufficient  size  The  ignition  question  is 
primarily  a  tbermodynamic  one;  the  likely  scenario  for  igniuon  based  on  a  dielectric  breakdown  event  is  as  follows 
In  each  local  breakdown,  some  amount  of  e.'iergy  is  dissipated  into  the  propellant  through  Joule  heaang  as  charge  is 
transfeaed  between  the  alumnum  pamcles  and  some  energy  is  also  dissipated  through  radiation  from  the  discharge 
arc  itself.  In  order  fev  ignioon  to  occur,  this  energy  must  be  dep^'cd  m  such  a  way  as  to  raise  the  temperature  of 
the  material  surrounding  the  arc  to  a  high  enough  valite  for  a  long  enough  time  for  a  self  sustaining  chemical 
reactioa  to  occur  m  the  propellant  How  high  this  temperature  must  be  and  how  long  it  must  be  maintained  is  a 
fiiDctioa  of  die  theimodynaimc/cheiiucal  properties  of  the  materials  mvolved  and  the  geometry  of  the  propellant  at 
the  microstnictuial  level.  These  are  issues  that  will  not  be  addressed  hue,  but  one  key  parameter  that  is  needed  to 
b^in  any  analysis  of  the  ignition  questitm  is  an  estimate  of  how  much  energy  is  dissipated  in  a  local  breakdown  or 
sequence  of  local  breakdowns.  Such  an  estimate  can  be  made  fiom  the  breakdown  simulations  described  m  Reference 
[1]  and  win  be  the  subject  of  this  section 

At  each  st^  in  the  sequence  of  local  bre^downs,  two  unportani  quantities  are  known  that  have  direct 
refevanoe  to  die  amount  of  enugy  released  dunng  such  an  event  Fust,  the  voltage  difference  between  the  metallic 
particles  is  known  just  befme  the  breakdown  event  frrnn  the  electrostatic  solution  for  the  potentials  In  addition,  the 
net  charge  on  the  particles  is  known  both  before  and  after  the  local  breakdown  what  the  new  electrostatic  solution  is 
foaod.  if  the  particles  were  not  part  of  any  imakdown  path  before  die  local  breakdown  in  question,  then  the  net 
charge  bef<ne  the  breakdown  was,  of  course,  zero  From  these  parameters,  the  energy  can  be  estimated 
stonghtforwardly  by  takmg  die  net  charge  cransfened  across  the  gap  between  the  particles  and  muluplying  that  charge 
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by  the  voltage  drop  that  ^sted  across  that  gap  before  the  breakdown  This  is  (xily  an  estimate,  of  course,  the  exact 
value  for  the  energy  released  is  a  function  of  the  exact  panicle  geometry  and  local  fields  This  is  an  order  of 
magnitude  estimate  only,  but  one  which  may  be  useful  as  a  startmg  point  for  more  complicated  and  detailed  ignition 
modeling. 

The  energy  released  m  each  local  breakdown  for  the  3D  simulauons  of  135  spheres  was  estimated  from  the 
simulation  data  as  described  above  The  maui  conclusum  is  that  the  energy  released  durmg  a  local  breakdown  between 
particles  is  on  the  order  of  10  to  100  picojoules  llus  estimate  comes  from  assummg  that  the  particles  have  a  radius 
of  microns  and  that  the  electric  field  between  the  particles  just  before  breakdown  is  nght  at  the  breakdown  field  of 
the  binder  which  is  taken  to  be  20  MV/m.  The  simulauon  data  can  be  scaled  to  give  values  for  the  energy  based  on 
any  particle  size  and  any  field  between  the  particles.  If  the  pamclcs  are  larger  that  20  microns  m  diameter,  then  the 
energy  released  is  larger  because  the  charge  transfer  is  larger  Also  if  the  field  between  the  particles  is  larger  than  the 
breakdown  field  of  the  binder,  this  also  mcreases  the  energy  released.  The  basic  relationship  which  scales  the  energy 
is 


U=CE2a* 

where  U  is  the  eneigy,  E  is  the  field  between  the  pamcles  the  apphed  field),  a  is  the  particle  radius  which  also 
determines  the  system  size  for  a  given  volume  fracoon  and  number  of  pamcles,  and  C  is  a  constant  of 
propomonaliiy  Hence  for  any  field  or  parade  size,  the  energy  can  be  estimated  fiom  the  simulation  data. 

Another  conclusion  suggested  by  the  simula&tm  data  is  that  more  energy  is  released  in  those  breakdowns 
that  occur  later  ui  the  breakdown  process  than  the  uuual  breakdowns.  Agam  this  assumes  Uiat  the  field  between  the 
particles  undergoing  breakdown  is  held  at  a  fixed  value  of  MV/m  regardless  of  when  they  break  down  This  result 
occurs  because  more  energy  is  available  in  a  region  where  the  field  between  two  particles  is  fixed  and  the  distance 
between  them  is  increased.  Recall  that  the  distance  betv^en  particles  that  undergo  breakdown  generally  increases  as 
the  toakdown  {voccss  coatuiues.  With  the  assumpuon  of  a  20  MV/m  field  between  parades  and  20  micron  diameter 
parades,  the  energy  rdeased  during  breakdown  in  samples  at  the  lowest  volume  factions  of  20%  ranges  from  5  pJ 
in  die  early  breakdowns  to  more  than  50  pJ  fOT  breakdowns  occurring  at  later  stages  m  the  b^akdown  process 

For  volume  fractions  of  50%.  the  highest  studied  fca  3D  systems,  the  energy  released  ranges  from  around  5 
pJ  fev  initial  breakdowns  to  only  10  pJ  later  in  the  breakdown  sequence  This  is  again  easily  interpreted:  at  higher 
volume  fractions,  the  avenge  parade  separation  distance  decreases  so  that  the  distances  involved  in  breakdowns  la^ 
m  the  sequence  are  still  quite  small  whereas  the  distances  buween  parades  m  the  mitial  breakdowns  are  quite  small 
m  eidier  case  A  final  comment  is  that  the  energy  stored  between  two  spheres  placed  m  a  known,  umf^m  field  can 
be  calculated  exactly  after  findmg  the  capaatance  between  the  spheres  In  the  case  of  a  20  MV/m  field,  the  energies 
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ace  on  the  onkrof  ins  of  pio^oi^  for  separaoon  distances  like  those  encountered  tn  the  breakdown  simulations 
This  is  not  ma^ecied  irrf  the  eoeigy  stored  tn  this  two  sphere  system  should  be  the  same  order  of  magmtude  as 
the  eaeigy  stored  betwcca  two  i^ibeies  in  a  random  ccKiflguraaon. 

In  niiinarj.  the  of  ignition  processes  in  solid  rocket  propellants  is  a  compbcated  problem 

involvix^  die  aaweoiop  of  stored  ekctcical  energy  to  heat  and  the  subsequent  ihermodynamic  transfer  of  this  beat  to 
the  materiai  nt^  or  inay  not  ignition.  Any  ignmon  modeling  therefore  dqiends  on  knowing  the 
magnitude  energy  rtieased  is  a  breakdown  process.  This  energy  is  on  the  order  of  10  to  100  picojoules  for  local 
hceskdowns  between  p— of  20  micruis  in  diameter,  dependmg  on  the  distance  between  the  pamclcs  which  is 
large*  at  fm^w  voScQse  farticmA  and  sm^iUer  at  targe  volume  fractions.  Hence,  although  propellants  with  a  higher 
volume  firaccim  of  break  down  easier,  more  energy  is  released  in  the  brealedown  of  propellants  with  a 

lower  vohsne  tacdon  of  ahuahMUL  Whether  this  amount  of  energy  is  large  enough  and  can  be  deposited  rapidly 
enoa^  (o  an  ugaidop  ^socess  to  b^m  is  a  question  m  need  of  further  investigation.  In  pan  it  is  covered  m 


CHAPTER  4 

IGNITION  OF  SOUD  PROPELLANTS 


4.1  Introduction 

Analyses  of  ESD-caused  ignition  (as  opposed  to  breakdown!  of  solid  propellants  must  include  an 
invesugation  of: 

1.  Elecincal  energy  available  to  Stan  igniiitHi  ui  typical  system-level  geometries, 

2.  Electneal  energy  deposition  in  typical  arcs  as  a  funcuon  of  sample  size  and  characteristics  such  as 
lesisiaoce, 

3.  Electcical  energy,  power,  and  time  requiremenis  for  arcs  to  cause  igmtion, 

4.  Resistive  loss  mechanisms  in  arcs,  and, 

5.  Physical  properties  of  solid  propellanis  (e.g„  heat  capacity)  that  relate  to  igmuon 

The  above  considerations  can  be  broken  up  into  two  levels  of  invesugation  a  system  or  macroscopic  level, 
and  a  local  or  microscopic  level  In  earlier  secuons,  we  discussed  the  elecuostaucs  of  system  level  ESD  charging 
Our  discussion  of  the  local  or  microscopic  level  is  contauied  in  Secuons  42  and  43,  and  will  be  broken  up  into  the 
arc  decuodynamics  and  the  igmuon  physics.  Sectton  42  addresses  the  relauonship  between  energy  and  propemes  of 
the  arc,  and  Section  43  discusses  the  basic  laws  relabng  to  propellant  igmuon,  and  the  esuapolanon  of  current 
igmuon  threshold  data  to  the  short  ame  scale  of  micro  discharges. 

We  note  diat  the  ESD  commumty  often  uses  the  words  vulnerabihty  and  suscepubihty  to  describe  aspects 
of  the  breakdown  and  ignition  problem.  In  our  termuiology,  the  available  energy  on  a  system  level  is  related  to 
vulnerability,  and  the  energy  in  the  micro  discharge  (between  particles  m  metal  loaded  propellants,  for  example)  is 
related  to  susceptibility. 

A  primary  objective  of  this  chqxer  is  to  provide  a  set  of  nomograms  to  quanufy  the  relauonships  between 
the  vanous  imponam  tgniuon  parameters.  There  are  many  parameters  to  keep  m  mmd,  includmg  power,  power 
density,  energy,  energy  density,  voltage,  charge,  current,  resistance,  capacitance,  ume,  and  charge  density,  for 
example.  There  is  general]/  a  physical  property  of  an  object  that  is  a  parameter  for  curves  on  the  nomogram 
Examples  of  such  pmameters  are  arc  length,  capaatance,  spark  g^  distance,  source  charge  density,  and  resistance  It 
is  hoped  that  drese  nomograms  will  help  the  reader  u>  appreciate  quanbtauvely  the  many  relauonships  that  exist  at 
the  system  and  local  levels.  Additional  informauon  that  may  be  helpful  m  understanding  the  concepts  in  this  paper 
are  contained  m  Reference  [2, 13, 14] 


4.2  System  Level  Electromagnetic  Principles 

At  the  system  level  we  will  show  the  relation  between  the  mtemal  voltages  and  fields,  and  the  external 
sources  which  are  charge  densuies.  The  physical  properties  of  the  system  that  relate  these  quantities  are  capacitance 
and  system  level  dimension.  The  mathematical  relatums  between  the  variables  are* 


Q  = 

QsA 

(4.1a) 

V 

Q/C 

(4.1b) 

V 

Ed 

(4.1c) 

^  = 

1/2  CV2  (orl/2QV) 

(4  Id) 

where 

Q  s  total  charge  (coulombs) 

Qs  s  charge  dcnsiiy  per  area  (coulombs  per  square  meter) 

A  «  area  (square  meters) 

V  B  voltage  across  the  system  (volts) 

C  s  capadtaocc  of  the  system  holding  the  source  charges  (farads) 

£  s  magmtude  of  electric  field  (volts /m) 

d  s  system  level  gap  (distance  over  which  E  IS  developed)  (meters),  and 

14  »  available  energy  Qoules) 

Figure  4. 1  shows  these  equations  in  graphical  form  (nomograms),  usmg  log-log  scales  for  each  nomogram 
Fev  example.  Figure  4.1a  (the  n^t-hand  figure)  shows  £q  4.1a.  where  we  can  enter  a  known  area  A  (abscissa)  and 
an  assumed  charge  density  Qs,  to  obtam  a  specified  charge.  These  nomograms  allow  one  to  more  clearly  visualize 
the  rdatiOD^p  of  the  many  system  level  electromagnetic  parameters. 

In  the  example,  we  assume  an  area  of  two  square  meters  and  a  charge  density  of  Qs^lS  ticoulombs  per 
square  merer.  This  gives  a  total  charge  Q  of  30  ticoulombs.  Travelmg  left  from  Figure  4.1a  to  Figure  4.1b  gives  a 
voltage  if  the  total  equivalent  capantance  is  known.  In  'a  example  case,  an  assumed  value  of  1  nanofarad  would 
lead  re  a  voltage  of  30  kilovolts  With  this  voUa^,  average  fields  (Figure  4.1c)  and  energies  (Figure  4. Id)  can  also 
be  estabhdied.  In  the  example,  a  field  of  1  MV  with  a  ^p  of  3  cm  as  shown  m  Figure  4.1c  gives  an  eno'gy  of  0  45 
joules  with  1  nF 
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Figure  4.1  System  Level  Relationships  Between  Available  Charge,  Voltage,  Field  and 
Energy  (See  text  f<x  defining  relationships  and  explanation  of  example) 


All  equations  here  were  hneai  except  the  quadmuc  dependence  of  the  available  energy,  Ua.  on  voltage,  as 
seen  in  the  steeper  dope  for  Figure  4.1d.  Figure  4.1  can  be  used  to  establish  upper  hmits  for  the  vanous  quanuues 
that  may  later  be  dictatmg  an  arc  and  possible  igniuon. 


4.3  Microstructural  Electro  Magnetic  Principles 

We  now  look  mote  closely  at  the  possible  values  for  discharge  energy,  to  assure  ourselves  that  they  arc  not 
in  fact  greater  than  the  available  energy.  Figure  4J2.  shows  four  interconnected  nomograms  Altliough  each 
nomogram  can  be  used  to  find  any  variable  as  a  function  of  the  other  two.  we  underline  below  the  most  likely 
variable  to  be  sdved  for  The  four  nomograms  respecQveiy  are: 

a.  Lower  lefc  Resistance  vs  voltage  with  current  as  a  parameter 

b.  tapper  left:  Power  vs  voltage  with  current  as  a  parameter 

c.  Lower  nghc  Resistance  vs  Time  with  capacitance  as  a  parameter 

d.  Uppernght:  Power  vs  time  with  power  density  as  a  parameter 


Figure  4.2  Discharge  Parameters  (Volts,  Current,  Resistance,  Power,  Capacitance,  Time, 
Energy) 
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'nie  equatioRS  for  these  nomograms  are* 


I  =  V/R 

(4.2a) 

P  =  IV 

(4.2b) 

t  =  RC 

(4.2c) 

U«Pl 

(4A!) 

By  eotenng  a  voltage  in  Fifure  AJa  (which  can  be  obtained  from  later  graphs,  with  knowledge  of 
breakdown  field  and  separation),  and  an  assomed  resstance  (obtained  from  independem  studies  of  breakdown  physics), 
we  can  detemune  the  likely  cnrrem  between  die  chscbjcgmg  spheres,  osuig  Eq.  4.'la.  In  the  example  shown,  we 
assume  10  volts  and  1000  amperes,  giving  a  cunent  of  10  milliamperes  (ma).  In  fact,  we  know  that  both  the 
voltage  and  the  resistance  are  dme-varying,  so  the  cocrent  certainly  is  also.  The  graph  should  only  be  used  to  get  an 
/<«»ge  value. 

Traveling  upward  with  the  same  voltage  to  Figure  4JZb.  and  osms  the  currents  just  obuuned  (altemauvely, 
we  could  have  used  the  resistance  with  die  formula  P  «  we  can  obtain  the  power  (time  rate  of  change  of 
energy  dqxsited),  us3g  Eq.  42b.  A  pos^le  error  can  creep  in  here,  since  the  cunrents  are  reversed  fiom  iheir  order 
inFigare42a.  In  the  exrmple,  we  obtain  a  power  level  of  0.1  Watts. 

In  the  lower  n^t  q^iadrant.  we  have  a  means  of  detenuining  the  time  scale  of  the  discharge,  using  Eq  4.2c 
The  time,  t,  is  read  off  of  the  abscissa,  with  the  discharge  resistaDce,  R  entered  os  the  ordinate,  usmg  an  assumed 
interqibere  capacitance,  C.  In  the  exanple,  we  take  O  0.03  locoforads,  so  the  assumed  R  of  1(XX)  ohms  gives  a 
discharge  tune  of  3  £*11,  or  30  picoseconds. 

In  the  upper  nght  quadrant,  the  results  obtained  fiom  Figure  42b  and  42c  are  combined  to  allow  an 
indqiendmit  d^emuoation  of  the  energy  deposied  by  tbe  are,  using  Eq.  42(L  In  the  example,  we  find  about  3 
{Hcojoules.  The  energy  so  obtained  should  not  be  larger  than  tbe  available  eneigy  obtained  from  the  charge  transfer 

Ud=  jivdt  =v|l(it  =  0^  QV  =  0JCV2  (4^) 

To  this  point,  we  have  estabhshed  a  system  levd  vuhiciabihty  energy  and  a  discharge  level  suscepubihty 
energy  (this  section).  Although  a  tnealolown  is  assmned  to  have  occurred,  tt  is  not  obvious  that  an  igniuon  will 
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ensue.  This  is  the  subject  of  the  next  secuon,  where  we  first  concentrate  on  those  physical  changes  (other  than 
ignition)  tha  will  first  result  hrom  the  breakdown 

4.4  Ignition 

Our  available  experimental  informauon  on  the  power  or  energy  density  required  for  igniuon  comes  from 
laser  and  BSD  experiments  A  major  difference  in  the  form  of  energy  appbcation  is  that  the  vast  majority  of  the 
laser  energy  can  be  applied  right  at  the  surface  which  must  be  heated  up  to  attain  an  igmtion  Only  reflection  of  the 
laser  beam  causes  mefllciency;  this  is  scHneumes  ccmtroUed  through  daricening  addiDves  Electnc  sparks,  on  the 
odier  band,  typically  demand  a  large  voltage  drop  within  the  propellant  or  along  its  surface,  which  is  only 
perqiherally  useful  in  heating  the  surface.  This  large  energy  may  possibly  be  minimized  by  using  low  resistances, 
but  this  in  turn  leads  to  shorter  arcs  and  a  higher  power  density  (but  not  energy  density)  requirement  Current 
breakdown  expenmenul  research  will  eventually  lead  to  a  good  understanding  of  this  trade>off  with  different 
resistances.  In  an  qieradon^  (not  laboratory)  situation,  there  is  a  selection  process  which  will  lead  to  those  being 
the  most  sensiove  naturally  occumng  first 

Another  difference  is  that  the  arc  diameter  can  be  much  less  than  the  several  milhmeiers  used  in  laser 
experiments.  For  small  currents  and  higher  pressures,  diameters  art  esumated  to  be  no  more  than  a  few  microns 
(hence  appreciably  less  than  a  millimeter).  However,  with  an  arc  to  a  surface,  the  arc  also  is  known  to  move  around, 
80  we  mi^t  first  stait  with  a  1  square  millimeter  as  a  possible  area.  In  laboratory  simulauons,  the  length  of  the  arc 
may  only  be  about  1  mm,  so  the  area  of  the  majenty  of  the  energy  deposition  may  be  even  less  than  1  sq  mm 

One  other  large  difference  for  electnc  spark-cau^  igmuon  is  that  the  enure  arc  is  at  high  temperature  not 
only  the  surface.  As  the  binder  and  ammomum  perchlorate  decomposes,  the  reacuon  products  will  be  ionized  and  (it 
is  believed)  can  thereby  more  readily  react  in  or  near  the  arc.  Smee  these  reacuons  are  exothermic  (they  release 
energy),  it  seems  possible  that  the  electnc  arc  may  be  more  efficient  than  the  lasers  m  causing  igniuon  Future 
research  will  possibly  'onfirm  this,  when  the  full  energy  balance  is  understood. 

Breakdown  condiuons  are  sometimes  cited  as  bemg  equivalent  to  ignition  in  terms  of  defining  a  senous 
safe^  hazard.  Indeed,  senous  concerns  have  been  rused  by  the  observation  that  igmuon  seems  to  have  been  well 
underway  (before  stepped  by  propellant  cracking)  without  a  full  discharge  across  test  samples  [Ref.  28]  There  is 
little  expenmemal  data  on  the  amount  of  energy  that  is  associated  with  these  anomalous  conoiuons,  but  the  required 
external  energy  would  appear  to  be  on  the  order  of  picojoules  However,  expenmemal  tesung  also  shows  that  there 
are  substantial  differences  -  certainly  m  most  experiments,  few  breakdowns  do  lead  to  igmuon.  However,  we  should 
note  that  in  high  pressure  test  configurauons,  other  mvesugators  do  find  that  most  discharges  will  lead  to  igmuon. 
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4.5  Micro-Discharge  Damage 

Here,  we  explore  obtaining  a  hazardous  situMion  with  micnwiischarges  (any  arcs  that  occur  between  small 
qiheies)  of  very  small  energy  content  Although  the  energy  content  is  low.  the  power  density  can  be  very  high, 
given  Uie  small  dimensKHis  of  an  arc  between  spheres.  We  will  explore  this  question  with  the  set  of  nomograms  in 
Figure  4.3.  The  primary  quanddes  of  interest  are  the  breakdown  field  and  discharge  energy  The  required  equauons 
am: 


ie=v/h 

(4.3a) 

U  =  ^cv2 

(4.3b) 

U^CydTvc 

(4Jc) 

Vc  =  p  h 

(43d) 

As2prh 

(4.3e) 

U<)  =  0/A 

(4  30 

where, 


£  s  electnc  field 

V  =  voltage 

h  B  height  of  discharge 

U  s  discharge  energy 

Ud  B  discharge  energy  density 

C  s  csfKicitance 

Cv  -  constant  volume  specific  heat 

dT  B  temperature  change 

Vc  B  arc  cylinder  volume. 

A  B  arc  cylinder  lateral  area 

Each  equation  is  {dotted  in  the  six  {^arts  of  Figure  4.3  Our  jinncipal  mterest  is  in  determining  the  radius 
of  the  arc.  assuming  that  the  breakdown  occurs  Solvmg  the  set  of  equauons  4.3,  we  find  • 
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r  =  E(Ch)-5/(6.28CvdT)5  (4.3g) 

However,  tbeeSects  of  ^  four  separate  parameters  (E.C,  hrCvdT)  are  hard  to  visualize,  and  the  nomogram  gives  a 
more  comptct  visnUzabon  of  the  reiationships.  Figure  43a  (lower  left  diagram)  allows  us  to  specify  an  assumed 
cyiinder  beis^  and  bieakdown  field.  In  our  first  example,  we  assume  that  an  average  field  of  about  1 0  MV  An  might 
be  aiffkaent  fdus  is  possibly  low),  giving  a  sphere  to  sphere  voltage  of  about  SO  voIls  This  is  an  appreciably 
laiper  voltage  tbao  leqatred  tor  breakdown  across  say  100  microns  of  aluminum  oxide  with  a  breakdown  field 
Sbcagfii  of  perhaps  1  GVAn.  where  only  about  10  volts  would  be  required. 

Ttavding  upwards  to  Figure  4  3b.  we  can  intersea  with  capacitive  lines  near  0.01  picofarads,  appropriate  to 
a  X  microfi  sphere  soe  (usmg  C»4rcea  for  a  single  sphere)  Even  with  the  closer  spacing  of  the  100  angstrom 
fanaiplftof  case  2  (aod  consequent  higher  capacitance),  we  iieduce  that  the  case  1  energy  is  higher  by  almost  a  factor 
of  10.  with  aa  eaeigy  of  123  picojoules 

We  JB&iBue  a  tomat  capaci^  of  10^  joules  per  cubic  meter  per  degree  Kelvin  and  a  desired  temperature  nse 
of  1000  degrees,  so  the  CydT  produa  is  10^*  This  value  of  thomal  capacity  needs  more  research  (that  is  currently 
onder  way),  bit  file  use  of  nomograms  allows  a  rapid  determinauon  of  :he  sensiuvity  of  this  parameter.  To  keep  the 
vais  managebbie,  we  conveit  to  picojoules  and  cubic  microns  (with  a  cubic  micron  denoted  as  cu  m  the  figure  and  1 
CttslO*^^  cibic  inetecs,inakmgtbeappropnate  diagonal  labeled  1000  pj/cu).  We  thus  predict  approximately  012 
aod  .0015  cidic  mkroits  (cu)  respecuvely  for  the  two  cas^ 

Tiivdiiig  downwards  to  part  d  of  Figure  4.3  to  find  a  cylinder  radius  (using  £q  4.3d).  we  find  about  0.03 
and  03  microas  respectively.  The  aspect  rauos  (h/r)  for  the  two  cases  are  very  differeni  -  about  167  and  0  02,  the 
fooner  a  irag  diia  cyfinder  and  the  latter  a  disk-shaped  object  The  temperatures  predicted  here  would  not  be 
aefaKved.  since  we  are  violaung  (he  assumpuon  of  a  relauvely  homogeneous  maienal.  with  the  large  alummum 
spheres  so  nearby  m  case  2. 

Pan  e  of  this  figure  Oower  nght)  cames  the  computauon  of  the  radius  (obtained  in  pan  d)  an  additional  step 
loobtainthclaieralsur6ccarcaofLhecylindncaIdjschargepath,usmgEq  A2&  The  lower  pan  of  this  figure  is  not 
vaSd  for  estunttmg  energy  transfer,  smee  the  height  of  the  cylmder  is  there  smaller  than  the  radius  However,  the 
Donogram  stiB  accnraiely  gives  the  lateral  surface  area  V/ith  this  qualifirauon,  the  "total"  area  for  dissipation  of 
energy  in  the  discharge  c  read  along  the  abscissa  of  Figure  4.3e  Note  the  closer  spacmg  of  the  r  diagonal  lines, 
since  tbe  lalesi  area  is  only  dependent  on  the  first  power  of  the  radius  (Eq  4.3e),  whereas  the  volume  of  Figure  4.3d 
is  d^ieodeniOB  the  second  power  In  the  example  showij,  a  lateral  area  of  3  square  microns  (abbreviated  su)  is 
fiwwn. 
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h,  spark  Unglh  (microns)  U,  Oap  Ensrgy  (plcojoulss) 


Traveling  upward  to  Figure  4.3f,  we  can  bnng  ihe  energy  obtained  in  Figure  4.3c  over  honzontally,  and 
fmd  at  the  intersection  the  energy  density,  using  £q.  4  3f.  in  this  example  the  energy  density  is  about  3  pic  jjoules 
per  square  micron.  This  value  of  energy  density  can  then  be  used  to  determme  the  likebhood  of  ignmon  using  the 
igmuon  chaiactensucs  of  the  next  secuon 

V,  Poienllai  Olflartnc*  (volts)  •P"*  volum*  (euble  mlctons)  Arot,  square  mierofts 


Figure  4.3  Characteristics  of  the  Discharge  Arc:  Volt.ige,  Field,  Height,  Capacitance, 
Volume,  Radius,  Area,  Energy  and  Energy  Density  (Surface  and  Volume) 

4.7  Ignition  CharacterisMcs 

Figure  4  4  introduces  the  subject  of  igniuon  chaiactensucs  Here  the  mterest  is  in  the  magnitude  of  the 
power  and  energy  densities  that  are  required  to  cause  igniuon  as  a  function  of  the  length  of  time  of  application  of  the 
power  or  energy  The  new  aspect  of  tnis  is  the  continuation  of  familiar  arguments  into  the  regime  of  sub 
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h.  ipitk  Inigm  (mictons)  U.  Gap  En.igy  (picoloul.i) 


nanosecond  pulses.  Here  the  anchor  value  is  taken  as  1  wait  per  square  mm  at  100  milliseconds  (based  on  data  in 
Hennance  [Ref.  2],  showing  approximately  25  calorics  per  c^2-sec,  and  pnvate  ccxnmumcations  with  workers  m  the 
field  [Ref.  IS],  givmg  about  the  same  value).  We  see  that  at  the  lowest  ume  scale  on  the  graph  (10'^^  seconds), 
iqi^oxtmately  ICO  kW/sq  mm  is  required  for  ignition  -  five  orders  of  magnitude  greater  than  for  a  time  of  0.1 
seconds.  However,  note  thtu  this  power  level  seems  more  manageable  when  considered  at  a  micron  level  -  it  is  only 
0.1  watts  per  square  nucron.  Figure  4.4a  also  shows  a  cross-hatched  region  of  departure  from  the  -2  slope,  due  to 
insufficient  pressure.  In  this  region,  an  mcrease  m  presaire  is  required  to  return  to  the  original  slope 

Figure  4.4b  shows  the  same  plot  as  it  is  displayed  by  mdividuals  interested  m  BSD  damage  of  electronic 
parts  Now  the  si(^  is  seen  to  have  a  value  of  -0.S.  Although  the  electronics  damage  specialist  also  talks  of  a 
slope  of  negative  .5.  that  person  is  talking  only  of  meUing,  and  of  course  not  of  igmuon  as  we  are  doing  here 
However,  the  similarity  m  slopes  and  presentation  is  of  course  because  both  are  explamed  as  thermal  phenomena 
The  change  in  slope  due  to  pressure  appears  different  in  Figure  4  4b,  but  is  of  course  due  to  the  same  phenomenon 
This  phenomenon  is  not  noted  in  the  eluctromc  parts  lirerature,  where  they  are  normally  looking  at  damage  deep 
inside  a  hgh  tesisdvi^  pcntion  of  a  mulu-layered  chip,  and  the  thermal  melting  is  unaffected  by  external  pressures 

The  presentation  reversal  of  Figure  4.4b  is  eq>ecially  useful  m  gomg  to  Figure  4  4c,  which  shows  the 
energy  m  the  same  igmaon  event  Most  importantly,  the  energy  required  for  igniDon  actually  decreases  with 
deocasing  tune  (or  increases  with  increasing  time).  As  the  product  of  power  and  time,  the  log  energy  vs  log  time 
plot  is  found  to  gave  a  slope  of  *.5  4  t  s  4.5  This  is  relauvely  easily  perceived  in  Figure  4  4c,  with  a  simple 
increase  m  the  slope  by  the  power  of  1,  due  to  the  muluphcanon  of  power  by  tune  to  obtain  energy  The  reducuon 
of  die  GO  area  due  to  insufficient  pressure  is  now  seen  as  a  horizontal  line  defining  a  region  m  which  reduction  of 
time  is  DO  longer  able  to  cause  igmtion,  as  it  is  above  and  to  the  nghl  of  that  line 

The  main  advantage  of  this  presentation  is  seen  in  Figure  4.4d,  which  is  the  cwnpanion  to  Figure  4  4a,  but 
with  a  sk^  (rf*  42  This  exponent  is  not  immediately  obvious,  nor  is  the  behavior  with  msufficient  pressure  The 
most  in^rtanl  feature  of  the  graph  are  the  very  low  values  of  energy  seen  with  short  ume  deposiuons  If  one  could 
get  all  die  way  to  point  2  on  this  plot,  then  a  depcsiuon  of  cmly  1  microjoule  per  square  nun  in  ^  seconds  is 
predicted  to  be  adequate  Expressed  in  turns  of  picojoules  and  square  microns,  this  is  the  same  number  1  picojoule 
per  square  micrcm. 

These  energy  and  temperature  values  seem  reasonable  in  terms  of  what  we  know  about  the  problem, 
although  they  are  much  smaller  than  expenmentally  reported  values  known  to  have  caused  an  igniuon  The  lowest 
Ignition  energy  we  know  of  is  on  the  order  of  mtUijoules  for  a  1  inch  sample  Microjoules  arc  reported  for  one  case 
when  extreme  pressure  was  also  present  We  think  that  a  lot  of  the  remaimng  difference  can  be  made  up  as  the 
energy  contained  in  muluple  breakdowns  mvolvmg  large  stnngs  of  parucles 
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The  remaining  variable  to  worry  about  is  the  energy  and  temperature  increase  required  for  the  Ammonium 
Perchlorate  particbs  to  participate  m  a  chemical  reaction  This  will  depend  greatly  on  the  size  of  the  AP  panicle, 
which  generally  is  much  larger  than  the  volume  calculated  above  for  the  discharge  volume  This  is  presumably  why 
sufHcient  energy  is  available  to  melt  (perhaps  vaptwize)  the  alumina  (and  a  small  amount  of  alummum),  but  not 
enough  to  cause  ignition. 


logP(w«n«K|mfn)  logtlmt(MCond9)  tog  U  oou»«s/*q  mm) 


Figure  4.4  Hypothetical  Propellant  Ignition  Characteristics 
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4.S  SUMMARY 


Hus  chapter  has  attempted  to  provide  a  brief  mlroducuon  into  present  understanding  about  the  ESD-caused 
Ignition  hazard  with  solid  propellants.  Although  the  understandmg  is  still  imperfect,  much  has  been  learned  since 
the  1985  Pershmg  n  accident  It  is  clear  that  the  hazard  can  be  ptimarily  attnbutcd  to  the  close  proximity  of 
aluminum  particles,  which  can  pteferennally  cause  a  breakdown  through  a  propellant  rather  than  around  it 

We  have  attempted  to  develop  a  graphical  mechanism  for  demonstraung  the  mter-relationships  between  all 
of  the  many  ESD  variables  •  both  for  ESD  and  for  igmuon  It  is  believed  that  these  nomograms  show  all  the 
vanoos  major  trends  more  clearly  than  would  tabulations  or  mulnple  graphs  with  parametric  vanauons.  The 
primary  cnbcal  variables  arc  the  diameters  and  spacing  of  pamcles 

Secondly,  we  have  explored  quanutauvely  the  possibUity  that  the  amount  of  energy  available  in  a  micro- 
discharge  might  be  sufficient  to  cause  an  ignition,  usmg  an  extiapolauon  of  experimental  results  obtamed  for  longer 
times  of  power  application.  It  seems  likely  that  the  energy  is  available  to  cause  material  melting  and  vaporization 
over  a  radius  of  about  one  thousand  angstroms  (0.1  micron),  when  the  gap  length  is  large.  More  work  is  required  to 
determine  the  likelihood  of  the  large  spaemgs  tequued  to  see  the  larger  energies  We  do  not  beheve  that  we  have 
proved  the  possibility  that  mtcnial  micro-discharges  are  adequate  to  cause  an  ignibon,  but  we  do  believe  that  we  have 
shown  it  is  possible. 

Primary  among  the  areas  requiring  more  research  are  more  accurate  measures  of  the  arc  resistance  More 
accurate  determination  of  pacameiers  such  as  the  duratton  of  the  micro-discharges  will  enable  the  resistance  to  be 
determined.  The  analysis  of  the  micro-discharge  is  found  u>  be  much  more  complicated  than  the  corresponding 
macro-analysis.  SigniScant  additional  reseaich  is  lequued. 
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CHAPTER  5 
SAFETY  OVERVIEW 


The  pievioiis  has  suntmanzed  the  final  anal>'tica]  work  at  EMA  dcahng  with  BSD  breakdown 

This  has  pnmanly  dealt  with  the  'genetic"  microscopic  analyses  that  we  have  performed  under  ARO  sponsorship 
However,  there  is  another  large  body  of  work,  both  at  EMA  and  elsewhere  that  has  dealt  with  specific  motor 
configurations  that  could  not  be  covered  in  this  study  and  this  report.  Although  much  can  be  learned  from  the 
general,  noa-spedfic  analyses,  the  specific  analyses  are  necessary  to  ensure  safety  for  specific  systems. 

Although  the  subject  of  safety  was  not  a  major  part  of  our  research,  it  has  an  overwhelming  importance  to 
the  industry  and  individuals  who  will  be  interested  in  this  report  Thus  we  feel  it  appropriate  to  close  this  report 
with  a  summary  of  what  we  know  about  safety  issues  identified  through  general  rniths,  which  are  summarized  m 
terms  of  the  four  mam  ESD  areas  described  m  Reference  2' 

Charge  Generation  Safety 

a.  Materials  should  be  selected  which  minimize  the  imtial  elecuostatic  generauon  Although  this  rule  is 

often  contradictory  with  other  requirements  (for  mmimum  fhcuon  for  example),  the  motor 
manulaictiiiing  commumty  has  been  able  to  find  replacement  materials  that  are  provmg  satisfactory. 

b.  Conductive  materials  are  generally  less  likely  to  generate  charge. 

c.  Proper  grounding  can  rrqitdly  remove  charge. 

Geometric  Field  Enhancement  Safety 

a.  Design  of  solid  propellant  motors  must  lake  account  of  the  ESD  hazard  m  the  future,  but  it  is  almost 
impossible  to  mtroduce  a  modiTicanon  into  an  exisnng  design  Proper  design  will  ensure  adequate 
qncmg  of  metal  pans  and  will  avoid  sharp  coiners  that  can  augment  the  mtenor  fields. 

b.  Metal  shieldmg  can  prevent  large  fields  through  the  Taiaday"  shieldmg  phenomenon.  This  is 
especially  unportant  for  lightmng  safety.  Enclosing  propellants  in  a  metallic  enclosure  is  the  oust 
means  of  ensunng  safely  The  use  of  caibon  fibeis  m  composite  cases  is  sn  excellent  altemauve 
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c.  >>4^ying  conductive  paint  to  nonmetallic  enclosures  is  of  high  importance  This  approach  has  been 
used  retroactively  with  several  motor  systems  and  has  been  under  crasiderauon  with  other  fielded 
systems.  Inspection  programs  should  be  in  place  to  ensure  that  these  paints  are  adequately 
ctmnecting  the  metallic  pans  of  the  motor  to  the  non-metallic 

d.  Grounds  must  be  in  place  to  tie  all  parts  of  the  motor  together.  Personnel  must  understand  the 
function  of  the  grounding  systems  through  training  programs.  In  some  circumstances,  personnel 
should  use  personal  wnst  straps  during  handlmg  operations. 

Microscopic  Field  Enhancement  Safety 

a.  The  size  of  the  alununum  and  other  particles  are  primarily  dictated  by  specific  energy  requirements 
and  considaatitms  othff  than  BSD.  Ifowever,  the  wide  range  m  breakdown  fields  among  different 
prt^ieUant  formulaticms  suggests  that  a  great  deal  of  BSD  safety  improvement  can  be  obiamed  if 
mme  attentioo  is  given  to  the  fonnulation.  The  immary  task  wQl  be  to  ke^  the  alummum  particles 
separated.  The  use  of  smaller  ammomum  perchlmate  particles  is  an  important  first  st^  m  this 
direction. 

b  The  conductivity  of  the  bmder  is  especially  important  in  establishmg  safety,  with  additives  being 
he^iful  to  increase  conductivicy. 

Ignition  Safety 

a.  Pressure  and  coafinunent  are  extremely  important  in  gomg  frcm  a  rdauvely  msigniiicant  l^eakdown 
to  the  more  cmastrc^ic  ignititHi.  Sups  must  be  taken,  e^iecially  durmg  the  manufacturing  jxocess, 
to  ensure  that  friction  is  low  and  that  mechanical  pressures  are  mmimized,  whenever  electnc  charges 
may  be  present 

b.  Since  discharge  time  is  so  unport^t  with  the  shorter  discharges  being  more  hazardous,  sufficient 
resistance  must  be  in  the  grounding  paths  as  to  avoid  a  very  rapid  discharge  Slow  discharge 
(microsecoQds  or  milliseconds)  ts  better  than  rapid  (picosecond)  discharge 
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